Molecular dynamics simulation of carbon nanotube growth under a tensile strain

We performed molecular dynamics simulations of carbon nanotube (CNT) to elucidate the growth process in the floating catalyst chemical vapor deposition method (FCCVD). FCCVD has two features: a nanometer-sized cementite (Fe3C) particle whose melting point is depressed because of the larger surface-to-volume ratio and tensile strain between the growing CNT and the catalyst. The simulations, including these effects, demonstrated that the number of 6-membered rings of the (6,4) chiral CNT constantly increased at a speed of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${1}\,{\textrm{mm}/\textrm{s}}$$\end{document}1mm/s at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${1273}\,{\textrm{K}}$$\end{document}1273K, whereas those of the armchair and zigzag CNTs were stopped in the simulations and only reached half of the numbers for chiral CNT. Both the temperature and CNT chirality significantly affected CNT growth under tensile strain.

The total MD simulation time was 500 ns , i.e., the CNT was pulled up for 0.5 nm , which is correspond to the length of the two 6-membered rings.We considered CNTs with diameters of 0.7 nm and the chirality of the armchair (5,5), chiral (6,4), and zigzag (9,0) CNTs to compare the differences in the growth processes depending on the CNT structure.All calculations were performed using the LAMMPS package 33 .There are many kinds of force fields for the mixture of iron and carbon systems, but almost all of them are not suitable for simulating the cementite structure, which is considered as the basic structure in the high-temperature FCCVD process.The previous study reported that only the embedded atom method (EAM) by Lau et al. 34 and Ruda et al. 35 and the short-range Tersoff-Brenner-type analytical bond order potential (ABOP) by Henriksson et al. 29 are the reliable potentials to describe the properties of cementite 36 .Among them, Henriksson's ABOP can reproduce the 6-membered ring structure of carbon allotropes 37 because it uses Brenner's reactive bond-order (REBO) potential 38 for carbon-carbon interaction, whose accuracy is well tested by comparison with density functional theory.It should be noted that the development of ReaxFF to simulate the mixture of iron and carbon system is continuing, but the simulation by ReaxFF potential takes much longer time than that by ABOP.Based on these references and our preliminary tests, we used the Henriksson's ABOP to simultaneously describe the interaction of atoms in the cementite and CNT in this work.

Results and discussions
Figure 2 shows the z-oriented trajectories of carbon and iron atoms in the cementite nanostructure at 1073 K , 1273 K , and 1473 K .The dynamics of the carbon atoms above the melting points ( 1273 K and 1473 K ) are well activated, and the carbon atoms are freely exchanged in their positions.The average velocity of the carbon in the direction normal to the surface was estimated using the following expression: The average velocities were 5.4 mm/s at 1073 K , 126 mm/s at 1273 K , and 143 mm/s at 1473 K , and all velocities were much higher than the reported velocities of CNT growth.Therefore, we can conclude that the bottleneck process of CNT growth is not the carbon supply from inside the cementite structure, but the C-C bond formation at the interface.The pulled-up speed for CNT growth in this study was 1 mm/s , which is relatively fast with respect to the experimentally reported speed, but quite a possible speed from the viewpoint of carbon sources.
Based on the above results, MD simulations of the growing CNT on the cementite surface were performed at a pull-up speed of 1 mm/s .Full trajectories during 500 ns simulations are available in Supplementary Video S1-S3. Figure 3a,c,e shows the final structures of the 500 ns MD simulations for armchair, chiral, and zigzag CNTs at 1073 K .This shows that the carbon supply from the cementite at 1073 K is insufficient, and the pulled- up CNTs are separated from the surface of cementite for all chiralities.Although the average velocity of carbon atoms in cementite was faster than the CNT pulled-up speed, only carbon atoms near the surface moved, and carbon diffusion from the bottom of cementite did not occur.Most of the atoms in the cementite maintained their initial positions, and the carbon atoms in the cementite were unable to catch up with the pull-up velocity of the CNT, resulting in the lifting of the cementite for the armchair CNT and the breaking of the connections between CNT and cementite for the chiral and zigzag CNTs.The trajectories of the number of 6-membered rings were nearly unchanged regardless of the CNT chirality, as shown in Fig. 3b,d,f, whereas those of the 5-membered rings increase with simulation time because dangling bonds of the CNT edge are formed and tend to form bonds with neighbor carbon to reduce the energy.In contrast to the low-temperature condition, the final structures at 1273 K indicate that the CNTs are connected to the cementite surface for all chiralities as shown in Fig. 4a,c,e.The number of 6-membered rings nearly monotonically increased from the initial value as the simulation time increased, and no stable defect rings, such as 5-membered or 7-membered rings were formed (Fig. 4b,d,f).Note that the increasing behavior of the 6-membered rings for the chiral CNT is different from that of the others at 1273 K .The number of 6-membered rings for the chiral CNT continued to increase throughout the simulation, whereas those for the armchair and zigzag CNTs were almost saturated at approximately half of the simulation time and only reached half of the numbers for chiral CNT.These results indicate that the chiral CNT can grow stably even at a speed of 1 mm/s at 1273 K , which is above the melting point of cementite.Under high-temperature conditions, although CNTs eventually elongate from initial length as shown in Fig. 5a,c,e, the CNT were absorbed onto the cementite surface, and the reconstruction of the CNT occurred in the initial stage of the simulations.Therefore, the number of 6-membered rings alternately decreases at the initial stage, increases below the initial value with the undulation of cementite, and finally increases over the initial value, as shown in Fig. 5b,d,f.The reconstruction of the defect rings during the CNT reconstruction process is shown in Fig. 5b.There is a possibility of defects owing to the high fluidity of cementite at 1473 K .Temperatures above the melting point are required for CNT growth at a speed of 1 mm/s , but too high temperature may cause defects, which change chirality or cause CNTs to separate from cementite surface.
To investigate the carbon diffusion in cementite, we analyzed the change in the z position of the carbon atoms in cementite over the simulation time.At 1073 K , only the carbon atoms at a depth of approximately 1 nm from the surface moved, and most of the carbon atoms in cementite almost maintained their initial position, as mentioned above.However, the carbon atoms moved over the entire movable area with a thickness of 2 nm at 1273 K and 1473 K because of the melting of cementite.Reflecting this difference in carbon diffusion, the distribution of carbon atoms after the 500 ns simulation also showed different behavior between 1073 K and above, as shown in Fig. 6.Regardless of CNT chirality, the number of carbon atoms decreased only near the surface at 1073 K , whereas the carbon distributions at 1273 K and 1473 K were flatter than those at 1073 K .At 1273 K and 1473 K , carbon atoms are supplied from the lower part of the cementite to the cementite surface, where carbon atoms have been reduced due to being used for CNT growth.Because more carbon atoms in cementite can contribute to CNT growth than those at 1073 K , CNT can grow at 1273 K and 1473 K at a CNT pull-up velocity of 1 mm/s ; in particular, chiral CNT can continue to grow.
Figure 7 shows the initial growth steps of chiral CNT.New 6-membered rings were obliquely constructed along the CNT edges.This indicates an apparent chiral growth, leading to a chiral (6,4) CNT, which is a neararmchair CNT and is easier to grow compared to pure armchair and zigzag CNTs, as reported by Artyukhov et al. 39 .The edges of the armchair and zigzag CNTs are parallel to the surface of the cementite structure; therefore, many carbon atoms should approach from the cementite to the CNT edge simultaneously for efficient CNT growth.However, the carbon atoms at the edge of the chiral CNT are not in the same plane, and thus a time lag for reaching the carbon atoms is permitted.Therefore, chiral CNT continue to grow stably at a speed of 1 mm/s at 1273 K , whereas the growth of the armchair and zigzag CNTs cannot completely catch up with the CNT pull-up velocity.

Conclusions
We performed MD simulations of the CNT growth process under the FCCVD method.Melting-point depression and tensile strain, which are features of the FCCVD method, were considered.Our simulation demonstrated that the carbon atoms in the cementite structure over its melting point diffuse well, and their speed is faster   www.nature.com/scientificreports/than 1 mm/s , which is the fastest speed of CNT growth in the FCCVD method.The CNT could not grow at 1073 K because of the insufficient supply of carbon, and could grow with many defect rings at 1473 K because iron atoms from cementite structure devoured a part of CNT structure.The number of 6-membered rings nearly monotonically increased in the chiral (6,4) CNT at 1273 K , whereas it was stagnant in the armchair (5,5) and zigzag (9,0) CNTs.

Figure 2 .Figure 3 .
Figure 2. Z position of the (a, b, c) carbon and (d, e, f) iron atoms in cementite at (a, d) 1073 K , (b, e) 1273 K , and (c, f) 1473 K during the last 50 ns of 100 ns equilibration.The colors from blue to red correspond to the z position at 50 ns.